In-situ accumulated stress measurements: application to strain balanced quantum 
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In this work we use the in-situ accumulated stress monitoring technique to evaluate the evolution 
of the stress during the strain balancing of InAs/GaAs quantum dots and quantum posts. The 
comparison of these results with simulations and other strain balanced criteria commonly used 
indicate that it is necessary to consider the kinematics of the process, not only the nominal values 
for the deposited materials. We find that the substrate temperature plays a major role on the 
compensation process and it is necessary to take it into account in order to achieve the optimum 
compensation conditions. The application of the technique to quantum posts has allowed us to 
fabricate nanostructures of exceptional length (120 nm). In situ accumulated measurements show 
that, even in shorter nanostrcutures, relaxation processes can be inhibited with the resulting increase 
in the material quality. 
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I. INTRODUCTION 

In-situ characterization techniques are one of the most 
powerful tools to control and monitor the kinematics of 
the epitaxial growth of heterostructures. They can give 
real time information about the evolution of the growth 
front and the formation of nanostructures, among others. 
The reflection high energy electron diffraction (RHEED) 
system is a standard piece of equipment in most molec- 
ular beam epitaxy (MBE) reactors, and there are many 
works about its usage to monitor changes in surface re- 
constructions^— the formation of quantum dots (QDs) 
and quantum wires (QWRs)^ optimum conditions to 
fabricate them and, even, about how to get the size 
and shape of those nanostructures from the study of the 
diffracted pattern. 5 In summary it is a mature technology 
with a long tradition in this field. 

Despite it was introduced back in the early 90's by 
Schell-Sorokin and Tromp, probably less known is the in- 
situ accumulated stress measurement (ASM) technique^ 
It basically consist on measuring the stress accumulated 
in a sample during the epitaxial growth by monitoring 
the changes in its curvature. The kind of information 
that can be extracted by this technique is very broad 
and ranges from the anisotropic strength of the surface 
reconstitutionsji the study of QDs, QWRs and quantum 
rings formation, 8,9 thermal expansion in heterostructures 
or the formation and evolution of dislocations and plastic 
deformation processes! 10 ' 11 

In its basic form, the bending of he sample is measured 
using the deflection of a laser beam on the sample surface. 
This method is easy to implement inside a MBE reactor 
since all the setup is outside the vacuum chamber. In this 
case, a lever-shaped sample is fix from one of its ends to a 
special sample holder. An aperture of enough size made 
on the center of the holder allows the lever to bend freely. 
This particular requirements of the holder and the sample 
itself is probably what prevents a general implantation of 
the technique in commercial MBE reactors. 



In general, two parallel laser beams hit the sample in 
a direction perpendicular to the surface, one on the fix 
end and the other on the free one. If substrate bends, 
we can measure the deflection of the beam that hits the 
free end compared to the other beam. Using two beams 
reduces the noise associated with mechanical vibrations 
and small temperature variations. The deflection can be 
recorded collecting the reflected laser beams with two 
segmented detectors^— 



II. WORKING PRINCIPLES 

Using the above mention geometry, the substrate cur- 
vature can then be calculated as (Stoneys equation): 



(d — do) cos a 
d H 2L 



(1) 



where d is the distance between the spots in the detec- 
tors, dO the initial distance between spots, du is the sep- 
aration between the laser beams, L the sample-detector 
distance, a is the incidence angle and A(l/i?) is the sub- 
strate curvature variation. If the deposited layer material 
has a lattice parameter larger than the substrate, then 
it suffers a compressive stress and the substrate bends, 
acquiring a convex curvature (Al/R > 0). On the con- 
trary, if the lattice parameter of the deposited layer is 
smaller than the substrate, the strain is tensile and the 
substrate becomes concave (Al/R > 0). 

Before the plastic limit, where the sample suffers from 
partial relief of accumulated stress through the forma- 
tion of dislocations, changes in the substrate curvature 
and the accumulated stress can be related by means of a 
modified version of the Stoney's equation, to include the 
biaxial character of the stress in the thin layers: 
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where R is the curvature radius, h is the thickness of the 
deposited layer, hs is the substrate thickness, a the stress 
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in th elayer and Ms=(l-vs) /Y s the biaxial modulus that 
relates the Young modulus (Yg) and the Poisson modu- 
lus (us). This equation is valid only under the following 
conditions: 

1. The thickness of the deposited layer and the sub- 
strate are much smaller than their lateral dimen- 
sions. 

2. The thickness of the deposited layer is much smaller 
than the substrate. 

3. The stress induced by the layer does not have a 
component in the direction normal to the sample 
surface. 

4. Substrate material is linearly elastic, homogeneous 
and isotropic. The deposited layer must also be 
isotropic. 

5. Edge effects are negligible and physical properties 
are homogeneous in planes perpendicular to the in- 
terface. 

6. The strain and shear deformations are negligible, 
in such a way that layer and substrate are within 
the elastic limit at all times. 

7. Substrate has no constraints to bend in neither of 
the two directions. This condition is no fully satis- 
fied in the described experimental setup. As we use 
lever shaped substrates with one ends fixed to the 
holder, we constrain the bending along the short 
side. If the lever satisfies b > 3a, with a and b 
the dimensions of the long and short sides, respec- 
tively, then the deformation in the transverse direc- 
tion will not influence the bending along the long 
axis and this condition can be fulfil. 

On the other hand, the crystal structure of the materi- 
als used in this work do not allow to fulfil the condition of 
isotropy. Eq.[5]is not valid and it must be adjusted to the 
experimental conditions, taking into account the crystal 
orientation of the interface and the elastic constants of 
the material in the direction along which the curvature 
is measure. In this way, the biaxial modulus becomes: 
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where c^ are the substrate elastic constants. Rearrang- 
ing Eq. [5] for uniaxial stress along the [110] and [1-10] 
directions: 
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The value of C44 is approximately Mg/2 so the error in- 
troduced by using Eq. [5] instead ofH]is less than 2%. 

Until now, we have considered the accumulated stress 
introduced by a layer in an static situation. However 



during a MBE growth the stress variation might be due 
to changes in the deposited layer thickness, changes in its 
stress or even the surface reconstruction. For this reason, 
if the thickness of a layer changes dh in a time t+dt, using 
a differential form of Eq. [2j 
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The right hand side of this equation have three terms. 
The first one describes changes in stress associated to 
an increase of thickness h in a time interval [t, t + dt]. 
The second term accounts for relaxation processes in the 
already deposited layer at the time t. Finally, the third 
term is related with changes in the surface stress. We 
can define the accumulated stress at the time t as: 



dh 



da 



Ea[h(t)]=a( Z = h,t)- + Iq ^ 



Substituting Eq. [5] in [5] we obtain: 
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As it can be seen, the magnitude measure in this ex- 
periments is the sum of the accumulated stress and the 
stress associated to changes in the surface reconstruction. 
Finally, combining Eq. [T] in its differential form (taking 
cos a = 1), and Eq. [7] we get: 

Z*[h(t)] + [Ars] = M f s[d{t }- d0] (8) 
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III. IMPLEMENTATION AND 
CHARACTERISTICS OF THE TECHNIQUE 

An important improvement of this technique, as it 
is implemented in the Institute de Microelectronica de 
Madrid (IMM), is the use of a large area CCD camera 
to record the two beams simultaneously. This method 
has several advantages over the segmented detectors. On 
the one hand, the optical alignment is considerably eas- 
ier since there is only one detector to be put in place to 
record both beams. On the other hand, it has larger dy- 
namical range, as the reflected spots are recorded at all 
times regardless of their separation and exact positions 
(within a reasonable range). Finally, it has comparable 
resolution to the segmented detectors method without 
the need of low noise amplifiers or other extra equipment. 

Figure ?? shows a detailed schema of the AS measure- 
ment system available at the IMM. The laser source (608 
nm) produces an intense beam that hits a beam splitter, 
leading to two perfectly parallel beams of similar inten- 
sity. The beams cross the optical window of the MBE 
reactor and reach the sample perpendicularly to its sur- 
face. One of the beam, hereafter reference beam (RB), 
hits the fix end of the sample so its reflection is not af- 
fected by the growth process. The other beam, hereafter 
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signal beam (SB), hits the free end of the sample and 
its reflection will be affected by the bending of the sam- 
ple and hence by the stress accumulated during growth. 
The measurement of the reflected beams is performed in 
a backscattering geometry, minimizing the error intro- 
duced by the approximation made in Eq. [8] 

The RB reaches directly the CCD whereas the SB 
crosses a prism to change its trajectory and send it to 
the camera. This prism is of capital importance in the 
setup and allows the usage of a CCD camera instead of 
the segmented detectors. Even in the case of the beams 
been reflected perfectly parallel, the distances between 
the spots would be of around 1 cm. In a more realistic 
case, where the beams diverge due to the deflection of 
the lever, the separation at a reasonable distance from 
the sample surface (~1 m in our case) could be of sev- 
eral cm, too large for most CCDs. Since the accumulated 
stress measurements depends only in the distance differ- 
ence between the spots and not on their absolute value, 
this approach does not have any effect in the results. 

Wc use a SpotOn CCD camera of Duma Optronics Ltd. 
and their acquisition software to get a beam positioning 
with sub-micron resolution. The distance between spots 
is sent to a custom software that converts it into accu- 
mulated stress, in real time, by means of Eq. [5] This 
software also records the opening and closing of the ef- 
fusion cell shutters, giving an exact match between the 
accumulated stress evolution and the materials growth. 
With this information, and assuming a typical distance 
of L = 975 mm between sample and detector, dn = 8 mm 
as the initial separation of the laser beams, and using the 
parameters characteristic of our substrates (GaAs, thick- 
ness hs—WO Mm, A/5 = 124 GPa), we obtain a maximum 
resolution of 0.02 N/m. This high resolution is normally 
not attainable due to vibrations and noise in the envi- 
ronment. Mechanical vacuum pumps, either from the 
MBE reactor or from nearby equipments, have the most 
detrimental effect and must be disconnected in order to 
perform high quality measurements. The real resolution 
in our system is normally between 0.05 and 0.1 N/m. 



IV. EXPERIMENTAL RESULTS 

We have used this technique to characterize the growth 
of strain balanced quantum dots and quantum posts 
(QPs). This kind of nanostructures have a three dimen- 
sional shape and, thus, the common equations used to 
calculate the optimum strain balanced condition can not 
be used, as the strain is inhomogeneous and the layer 
thicknesses are not well defined.— ~— In this work we 
use the in-situ accumulated stress measurements, as de- 
scribed above, to obtain the real stress that introduces 
the QDs and the most appropriate GaAsP thickness and 
composition that exactly compensates that stress. 

All samples have been growth using solid source MBE 
on GaAs (001) substrates 100 /Ltm thick. InAs and 
GaAs/GaAsP growth rates are 0.02 and 0.5 ML/s, re- 
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FIG. 1. Stress accumulated during the formation of QDs after 
2 ML of InAs. The lower part of the figure shows the materials 
that are growing at each time. 



spectively. As beam equivalent pressure (BEP) is kept 
at 1.5xl0~6 mbar at all times. 



A. Strain balanced InAs quantum dots 

In Fig. [1] we show the evolution of the total accumu- 
lated stress as we grow an InAs QD layer for different 
substrate temperatures. As it was found by Silveira et 
al, four regions can be distinguished: 8 Region I: InAs be- 
gins to grow layer by layer, increasing the compressive 
stress linearly (except for a transition region at the be- 
ginning); Region II: just at the critical thickness, surface 
relaxes and QDs nucleate. The remaining deposited In 
keep floating on the surface or incorporates to the exist- 
ing islands but without increasing the stress; Region III: 
during capping this remaining In incorporates, suddenly 
increasing the accumulated stress; And Region IV: when 
In is exhausted, GaAs grows without any further change 
in the stress. As it can be seen, the maximum accumu- 
lated stress depends strongly on the substrate temper- 
ature and also on the total amount of In deposited, as 
shown in Fig. [UJl) and (b) (filled symbols). This kind 
of dependence is disregarded in the strain balance cri- 
teria used in QWs. The accumulated stress introduced 
by a flat, strained layer, assumed homogeneous, can be 
approximated by: 

Ecr L = M L e L t L (9) 

where M L , e L = {a su bs - aL)/asubs and t L are the layer 
biaxial modulus, the lattice mismatch between the layer 
and the substrate and the layer thickness respectively. 
Fig. [2] shows also the results of this equation applied to 
the nominal InAs thickness used in each case (open sym- 
bols). It can be seen that using the above equation to 
calculate the stress introduced by the QDs and, hence, 
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FIG. 3. Evolution of the accumulated stress as a function of 
the In composition for a given total In amount. Right scale 
shows the corresponding layer thickness. 



the strain balanced condition, leads to sub-estimate the 
accumulated stress in all cases. 

The reason for this discrepancy lays on the assumption 
that all the deposited In incorporates in the form of InAs. 
It is well known that, during the QDs capping, there 
is a large Ga-In intermixing, leading to quantum dots, 
wetting layer and capping made of InGaAs of varying 
composition. If the total In incorporated into the sample 
is to be constant, an InGaAs layer with a dilute alloy 
introduces more stress that a pure InAs layer. Fig. [3] 
shows the dependence of the accumulated stress on the 
In content of the layer (x), keeping the restriction of equal 
overall In content: 

tL x x = A = constant =>• £<tl = ML(x)ei,{x)aL,(x) — 

x 

(10) 

where t^ is the layer thickness in ML. M^(x) and aj,(x) 
are obtained as a linear interpolation of the GaAs and 
InAs parameters. As an example, if the In contained in 
a pure InAs monolayer is spread in two monolayers, giv- 
ing a Ino.5Gao.5As alloy, the resulting accumulated stress 
changes from -1.7 N/m to -2.1 N/m. 



FIG. 4. Accumulated stress as a function of the P BEP. Next 
to each curve is the P content of the alloy as estimated from 
Eq.! 



A similar analysis can be performed for the growth of 
the GaAsP compensating layer. In this case, the variable 
parameter is the P BEP, which gives the GaAsP composi- 
tion. Fig.|3]shows the evolution of the accumulated stress 
as a function of time. As it is expected, the accumulated 
stress in this case is tensile, owing that the lattice pa- 
rameter of GaAsP is smaller than that of GaAs. The 
composition of the layer can be estimated using equation 
[5J although it is not really needed for the calculation of 
the optimum strain balance condition. 

Knowing the compressive stress introduced by the QDs 
and the tensile stress compensated by the GaAsP strain 
balanced layer (SBL) as a function of its composition, we 
designed two strain balanced QDs stacks (A and B) aim- 
ing to a 100% of strain compensation. In both cases we 
use 2 ML of InAs for the QDs and a total spacer between 
layers of 15 nm. Substrate temperature and As BEP are 
kept constant at 510 °C and 1.5x10^6 mbar respectively 
during the growth of the stacks. The substrate is GaAs 
(001) cantilever shaped (4x20 mm) with a thickness of 
100 /im. The only difference between the samples is the 
compensating layer thickness and composition. In sam- 
ple A we use a 13 nm thick SBL with 4.3% of P after 1 
nm of GaAs capping, whereas in sample B we use a SBL 
5 nm thick and 18% of P after 8 nm of GaAs capping. 
The evolution of stress during the growth of both stacks 
can be seen in Fig. [SJa) and (b). 

Several things can be observed in these curves. Firstly, 
the average strain has been successfully balanced in both 
cases. Assuming that each QD layer introduces a stress 
of 5 N/m, we obtain an average strain compensation of 
95% for the sample A and 105% for sample B. Secondly, 
in the sample A, the oscillations corresponding to the 
accumulation/compensation sequence are damped. This 
is due to the InAs and GaAsP intermixing and was ex- 
pected given the small GaAs capping on top of QDs. The 
formation of quaternary InGaAsP compounds introduces 
a stress (compressive or tensile) smaller than the corre- 
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FIG. 5. Accumulated stress of samples A (A) and B (b). On 
the left there is an schema of the samples layer structure. 



sponding InAs or GaAsP alloys separately. Although the 
strain is balanced on average, the stoichiometry of the 
stack is uncontrolled. This intermixing has an impact on 
the optical properties of the QDs and must be taken into 
account when placing the barrier too close to the QDs. 

It should be notice that the results presented here are 
only an example of perfectly balanced stacks using the 
in-situ accumulated stress measurements. Other combi- 
nations of SBL composition and thickness are possible, 
such as using pure GaP layers, GalnP or dilute nitrides, 
having optical or electrical properties more suitable for a 
particular application. 



B. Strain balanced quantum posts 

Quantum posts (QPs) are assembled by epitaxial 
growth of closely spaced quantum dot layers, modulat- 
ing the composition of a semiconductor alloy, typically 
InGaAs. Contrary to normal self- assembled nanostruc- 
tures, the height of the QPs can be controlled by the 
number of periods of the superlattice grown on top of 
the seed QDs layer. The amount of In in this kind of 
nanostructures is very large compared to stacked QDs, 
with the result that the accumulated stress is enormous 
and there is a tendency to the formation of dislocations. 
The largest QPs reported are about 40 nm high (Pendi- 
ente de revisar). 

In this work we monitor the evolution of the accumu- 
lated stress in two QPs samples. The first one (sample C) 
uses a superlattice of 2.2 A of InAS and 8.5 A of GaAS 
grown at 510 °C. In sample D, on the other hand, we 
substitute the GaAs by GaAsP with a nominal 14% of P 
content. This approach has allowed us to fabricate ex- 
tremely large QPs of up to 120 nm with very interesting 
optical and electronic properties, as it has been reported 
elsewhere J£ 

Fig. [6] shows the resulting accumulated stress curves 
for both samples. Only the first 17 periods are mea- 
sured since the sample bending became to large to be 
recorded with the CCD camera. For each curve, there 
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FIG. 6. Accumulated stress of regular and strain balanced 
quantum posts (SB-QPs). Dashed lines are linear extrapola- 
tions of the accumulated stress in the first five periods in each 
case. 



is an approximated linear extrapolation of the stress ac- 
cumulated by the first five periods, marked with dashed 
lines. Grey areas in the background represents the peri- 
ods when the In effusion cell is open. The intermediate 
white regions represents a 10 s growth interruption under 
As flux plus the GaAs (GaAsP) growth. 

As it can be seen, sample C accumulates larger stress 
than sample D. The degree of compensation in the su- 
perlattice (disregarding the QDs seed) can be estimated 
from the linear extrapolations mentioned above and gives 
a value of 57%. The evolution of the accumulated stress 
in sample D is linear, behaviour that might be expected 
if all periods introduce the same amount of stress. The 
strain balanced oscillations are barely visible. As ob- 
served previously for QDs, this strong damping is directly 
related with the intermixing of the constituent materials. 
In sample C, it is remarkable the progressive bending of 
the accumulated stress that deviates from the linear ten- 
dency observed in sample D. Moreover, as it can be seen 
there is an apparent inversion of the accumulated stress 
during the In growth. In this periods, growing In reduces 
the accumulated stress, rather than increasing it. 

Both effects might be explained in terms of an ini- 
tial stage of relaxation processes in the superlattice. 
As shown by Ujue et al. during the growth of thick 
Ino.2Gao.8As layers on GaAs, prior to the formation of 
dislocations, there is an initial relaxation stage consist- 
ing on a ripening of the growth front along the [110] di- 
rection. From the point of view of accumulated stress, 
this effect produces a progressive deviation of the linear 
behaviour stated in Eq. [9j The onset of the relaxation 
depends on the growth rate, taking place earlier for slow 
growth rates. This is roughly the situation found in sam- 
ple C, where the average composition of In is also around 
20% and with an average growth rate of 0.07 ML/s. 

On the other hand, the reduction of the accumulated 
stress during In growth could be related also with a relax- 
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ation process but of a more local nature. The fabrication 
of QPs relays on the effective migration of In adatoms 
towards the top of buried QDs, where their elastic en- 
ergy is smaller. The accumulation of this In atoms could 
partly relieve the stress of the InAs beneath them by lo- 
cally increasing the lattice parameter of the structure. 
The growth of the GaAs capping suppresses this effect, 
increasing the stress by incorporating the InAs on the 
surface to the GaAs lattice structure. 

The comparison of sample C and D also indicates that 
using a strain balanced technique to grow QPs is not only 
necessary in the case of extremely large QPs, but that it 
could be desirable also in average size nanostructures, 
with more than 8 or 10 periods, to avoid the relaxation 
processes described above. 

V. CONCLUSIONS 

In this work we have shown the implementation of a 
very compact in-situ accumulated stress measurement 
setup based on the usage of a high resolution CCD 
camera to monitor the bending of the substrate during 
growth. The system outperforms previous designs in res- 



olution and simplicity. We have used this system to study 
the strain balanced process of QDs and QPs. We have 
found that it is possible to achieve perfect strain compen- 
sation in QDs stacks by calibrating separately the stress 
introduced by the QDs and the compensating layer. This 
process depends strongly on the substrate temperature 
and the incorporation of In atoms to the sample. Finally, 
we have shown a reduction of 57% in the stress accumu- 
lated during the growth of QPs by incorporating P to the 
matrix. The strain balanced technique is found to sup- 
press the relaxation processes that take place in the first 
stages of the grow of this nanostructures. In summary, 
this experiments show the capability of strain balance 
technique to improve the quality of quantum nanostruc- 
tures and the importance of kinematics in the optimiza- 
tion of the optimum strain balanced conditions. 
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